Dielectric anomaly in coupled rotor systems 
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The correlated dynamics of coupled quantum rotors carrying electric dipole moment is theoreti- 
cally investigated. The energy spectra of coupled rotors as a function of dipolar interaction energy 
is analytically solved. The calculated dielectric susceptibilities of the system show the peculiar 
temperature dependence different from that of isolated rotors. 
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I. INTRODUCTION 

With the advent of nanotechnologies, quantum rotors 
have attracted much attention in relevance to a funda- 
mental element of molecular scale machinery [1-3]. Ar- 
rays of surface mounted quantum rotors with electric 
dipole moments are of particular interest because dipole- 
dipole interactions can be controlled and even designed to 
yield specific behavior, such as ferroelectricity. Ordered 
two-dimensional arrays of dipole rotors yield either ferro- 
electric or antiferroelectric ground states, depending on 
the lattice type, while disordered arrays are predicted to 
form a glass phase [4, 5]. 

Besides technological problems, the microscopic dy- 
namics of quantum rotors have extensively studied from 
physical and chemical interest. The idea of quantum ro- 
tors is applicable to interstitial oxygen impurities in crys- 
talline germanium, where oxygen atoms are quantum- 
mechanically delocalized around the bond center posi- 
tion [6] . The rotational of oxygen impurities around the 
Ge-Ge axis has been experimentally observed by phonon 
spectroscopy [7]. While the rotation of oxygen impu- 
rities in Ge is weakly hindered by an azimuthal po- 
tential caused by the host lattice, several materials are 
known to show a free rotation of molecules. An exam- 
ple is ammonia groups in certain Hofmann clathrates 
M(NH 3 ) 2 M'(CN) 4 -G [8-10], usually abbreviated as M- 
M'-G, where M and M' are divalent metal ions and G is 
a guest molecule. Nearly free uniaxial quantum rotation 
of NH3 has been observed for the first time in Ni-Ni- 
(CeD 6 )2 by inelastic neutron scattering [8]. Recently, a 
surprising variation of the linewidth has been observed 
for Ni-Ni-(Ci2Hio)2 [H]> which has been interpreted by 
a novel line broadening mechanism based on rotor-rotor 
coupling [12]. It is also known that the (3 phase of solid 
methane [13] as well as methane hydrate [14] show al- 
most free rotation of CH4 molecule. The linewidths of 
methane in clathrates show inhomogeneous broadening 
owing to the dipolar coupling with water molecules [15]. 
It is therefore expected that new interesting phenomena 
will be found by investigating the influence of dipolar 
interaction between quantum rotors. 

In the present paper, we study the correlated dynamics 
of coupled quantum rotors carrying electric dipole mo- 
ments. We give the exact solution of eigenvalue problem 
of interacting rotors with arbitrary configurations. It is 



revealed that coupled rotors show a peculiar dielectric re- 
sponse at low temperatures, which can be interpreted by 
taking account of the selection rule of dipolar transition 
for coupled rotors. 



II. THE HAMILTONIAN 

Suppose two dipole rotors qr\ and qr?, separated by 
the vector R. The Hamiltonian for the system is given 
by H = Hk + Wd , where the kinetic term is 
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Here / is the moment of inertia for dipole rotors and s 
the dielectric constant, respectively. Figure 1 shows a 
configuration of two dipoles rotors under consideration. 
We assume that rotors do not feel any potential variation 
along a ring of radius r. In the Jacobi coordinate, the 
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FIG. 1: Schematic configuration of coupled rotors. Each ro- 
tors represented by ri and r-z rotates along a ring of a radius 
r, and carry dipole moment gr 1 ( 2 ). 
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9 2 )/2, Eqs. (1) and (4) yields 
H K = 
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FIG. 2: Contour plot of the interaction term Wd in (Oi, O2) 
plane. Two maxima (white regions) and two minima (dark 
regions) are realized at positions with differences AOi w it 
and AO2 ~ 7r. Parameter values are given in the text. 



vectors n , r 2 and i? are given by 
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A spatial profile of Wd as a function of (©1,02) is dis- 
played in Fig. 2 by a contour plot, in which the angles 
(ai,a 2 ,(3) are set as (7r/4, — ir/6, n/3). We should remark 
that two minima (dark regions) and two maxima (white 
ones) are located at the anti-parallel or parallel dipolar 
configuration, indicating that the dipoles prefer an anti- 
parallel configuration. The two minima of Wd(@i,@ 2 ) 
arise from the dipole interaction between two rotors, i.e., 
the dipole interaction plays a key role for creating barri- 
ers and two potential minima, which strongly affect the 
energy spectra and the dielectric response of the system. 

Provided that the spacing R is large enough compared 
with the radius r, the interaction term Wd can be ex- 
panded in terms of 1/R. The lowest-order term has the 
form of a dipolar interaction given by 
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The higher-order term wffl = Wd — W^ is of the order 
of 0(r 3 /i? 4 ), which can be negligible for the case R^> r. 
Actually we have confirmed that the calculated results 
presented in this paper change very little by taking into 
account the term Wp\ 

III. EIGENVALUES AND EIGENFUNCTIONS 

The Schrodinger equation for the Hamiltonian Hq = 

(3) 

Hr + W D has analytic solutions as shown below. Trans- 
forming variables to 9\ = (61 + 62)72 and 9 2 = (61 — 



The parameters a and 7 i (i = 1, 2) are functions of angles 
ai,a 2 and (3 defined in Fig. 1, whose explicit forms are 
given by 



c * = 7J\A+yf> 7* = ^ tan 1 ^ _ 
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with the definitions 



xi — sin/3(cosai — cosa^), 

x 2 — sin/3(cosai + cos«2), 

2/i = cos/3(l — cosai cosa 2 ) + 2sinai sina2, 

y 2 = cos/3(l + coscti COSC12) — 2sinai sina2- (8) 

Consequently, we can decompose the Schrodinger equa- 
tion i7o^o(#i, #2) = E <Fo(Si,S 2 ) into two independent 
Mathicu equations. Setting ^0(^1,^2) = </?i(#i) ¥'2(^2), 
we obtain 



+ —[c i E D cos2{6 i +~t i ) 



Ei]<pi = 0, [i = l,2] 
(9) 

where the quantities Ek = h 2 /(2I) and Ed = 
q 2 r 2 1 (4tteR 3 ) represent the kinetic and interaction en- 
ergy, respectively. The eigenvalue E of the initial 
Schrodinger equation is expressed as the sum of E = 
Ei + E 2 . Note that the periodic terms oc cos2(6>i + 7^) 
originate from two minima (or maxima) of the interaction 
term ^(61,62) shown in Fig. 2 [17]. 

Eigenfunctions of Eq.(9) are described by four 
types of the Mathieu functions, given by ce 2n (vi,6i), 
se 2n+ i(v l ,9 l ), cc 2n+1 (v i ,8 l ) and sc 2n+2 (v l ,8 l ) with the 
definitions i>, = CiEd/Ek and n = 0, 1,2- ••. Each of 
them belongs to a different eigenvalue and can be ex- 
pressed in terms of the Fourier-cosine expansion; for in- 
stance, 

00 

ce 2n (vi,6i) = 4™ } (^)cos2m(^ + 7i ). (10) 

m=0 

The coefficients {ifj'} are determined by a successive 
relation obtained by substituting Eq.(10) into Eq.(9). 

The amplitudes of {^m"*} rapidly decrease with increas- 
ing m, so that we can truncate the summation in Eq. (10) 
at m — 20 in actual calculations. 

Figure 3 plots the calculated spectra of eigenenergies 
E = Ei + E 2 as a function of Ed, where Ek is taken 
as an energy unit. The angles (a\,a 2 ,(3) are set to be 
(•7r/4, 0,0) as an example. We find, though some levels 
are degenerate when Ed = 0, they split off for finite 
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FIG. 3: The energy spectra of the paired-rotor as a function of 
Ed- Solid arrows indicate a part of allowed dipole transitions 
for the component p x in the case of (ai,a 2 ,/3) — (tt/4, 0, 0). 
The explanation on three labels (a)-(c) is given in text. 

Ed with a monotonous variation with increasing Ed- 
For high- .Ed limit, some levels become degenerate again. 
It indicates that the relative motion of paired-rotors is 
frozen out for Ed 3> Ek due to the strong Coulomb 
interaction. This behavior can be understood from the 
spatial profile of the interaction term Wd(@i, ©2) shown 
in Fig. 2. With increasing Ed, the depths of two minima 
of IFd(6i,92) grow, and larger barrier-heights hinder 
the quantum transition of a particle through the bar- 
rier. This gives rise to localized wavefunctions around 
these minima. Consequently, in the limit of Ed S> Ek, 
the amplitude of the eigenfunctions are strongly local- 
ized around two minima, and these two localized eigen- 
states are nearly degenerate. Even if the higher-order 

term wffl is taken into account, the energy spectra docs 
not change much, since it only slightly disturbs the sym- 
metry of the depths of two minima shown in Fig. 2. 
When varying the angles (ai,a2,(3), the curves in Fig. 
3 slightly shift to upwards and/or downwards except for 
the unchanged values of E at Ed = 0. 

IV. DIELECTRIC SUSCEPTIBILITIES 

Let us consider the dielectric response of dipole ro- 
tors coupled via dipolar interaction. The real part of the 
frequency-dependent dielectric susceptibility is expressed 
as 

X^,T) = |(^KI^)| 2 

Ej —Ei ( Ej \ 
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FIG. 4: The dielectric susceptibility x xx (T) for the zero- 
frequency limit u> — > as a function of the inverse temperature 
1/T. The strength of dipolar interaction is increased from 
top to bottom; i) Ed = (solid), ii) Ed = 0.01 (dashed), iii) 
Ed = 0.1 (dash-dotted), and iv) Ed = 1.0 (dotted) in units of 
Ek- Bumps at around £if/(fcsT) w 5.0 appear in the cases 
of ii) and iii) . Inset shows three components of r\^i (T) for the 
case of iii), whose definitions are given in text. 



where Z = ■ exp(— Ej/ksT) is the partition function, 
and \Ej) is the eigenvector belonging to the eigenvalue 
Ej. The quantity is the /x-component of the total 
dipole moment p = q(r\ + rz), which depend on the rel- 
ative orientation with respect to the external field. We 
should note that the susceptibility depends on the selec- 
tion rules for dipole transitions between different eigen- 
states. In Fig. 3, allowed dipole transitions for p x are 
indicated in part by solid arrows. Note that only a part 
of allowed transitions are shown in the figure, which are 
dominant for the dielectric susceptibility x(u, T) at tem- 
peratures T « Ek I ks ■ The rest of allowed dipolar tran- 
sitions do not contribute to the susceptibility given by 
Eq.(ll), because the energy difference \Ej — Ei\ is so 
large and/or the Boltzmann factor exp(— Ej/ksT) be- 
come much smaller than unity. The interpretation on 
three labels (a)-(c) shown in Fig. 3 will be given later. 

We have calculated the temperature dependence of the 
dielectric susceptibility x^fu, T) for various Ed- Figure 
4 shows the calculated results of dc susceptibility \xx (T) 
normalized by a factor q 2 r 2 /{eEk)- We have taken four 
values of Ed/Ek] the solid line (Ed = 0), the dashed one 
[Ed/Ek = 0.01), the dash-dotted one (E D /E K = 0.1), 
and the dotted one (Ed/Ek = 1.0). The angles are set 
to be (oti,ot2,l3) = (7r/4, 0,0) for all Ed- For the case of 
Ed = 0, the susceptibility monotonically increases with 
decreasing temperature, and becomes constant at lower 
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temperatures. The crossover temperature between the 
steady increase and the almost constant value in Fig. 
4 is determined by the minimum-energy difference of 
eigenstates at Ed = that are allowed for dipole tran- 
sition, namely, indicated as (a) in Fig. 3. For the case 
of Ed/Ek > 1, the strong Coulomb interaction prevents 
from the relative motion of rotors so that the magnitude 
of the susceptibility x(T) decreases with increasing Ed. 

It is noteworthy that, for relatively weak interaction 
Ed/Ek < 0.1, a bump is appeared in the suscepti- 
bility at about E K /(k B T) w 5.0. The kinetic energy 
Ek — h 2 /(2I) for actual rotating molecules is the order 
of 1 mcV [14], indicating that the characteristic temper- 
ature T* = Ex/ks x 0.2 corresponding to the bump is 
estimated as about 1 K. We made sure that the bump 
can be observed for any angles (a\, a 2 , 0) when Ed/Ek 
is less than 0.1. This anomaly stems from the correlated 
rotation of paired-rotors via the dipolar interaction, and 
can be interpreted by the argument on the selection rule 
for dipolar transition. 

To understand the origin of the bump, we decompose 
the total susceptibility Xxx(T) give in Eq.(ll) as 



2 K^bxl^i) 



(12) 



eZ Ej - Ei 



cxp 



E 3 
k B T 



exp 



Ei \ 
k B Tj 



(13) 



where Ylrj xfj\ is the summation over all possible com- 
binations of (j, I) under the condition I ^ j. Note the 
fact that only three components of i]jj(T) are responsible 
for the total susceptibility (12) around the characteristic 
temperature T*. We denote those components by r) a , 
and rf , which are characterized by the eigenfunction 
= (6 1 ,6 2 \E j ) and $ t = (6 1 ,6 2 \E l ) as follows; 



n 

*4 



*j = ce o (0i)ceo(0 2 ), *i = cei(0i)cei(0 2 ), (14) 
!Pj = ce o (0i)cei(0 2 ), *i = cei(0i)ceo(0 2 ), (15) 
rf; = sei(fli)se 2 (fl 2 ), #i = se 2 (0i)s ei (0 2 ). (16) 



The alphabets subscribed on 77 correspond to three dipo- 
lar transitions labeled by (a)-(c) shown in Fig. 3. For 
example, the solid arrow of (a) in Fig. 3 connects the 
eigenstates <Pj and <Pj defined in Eq.(14). 

For weak coupling Ed -C Ek, the solution of the 
Mathieu equation (9) is easily solved. In the lowest order 
of the perturbation theory, the eigenvalues Ei(i = 1,2) 
read in 

Ei = ^n 2 +aE D (n = 0, ±1, ±2, • • • ) (17) 



with a constant a. The solution (17) gives the eigenen- 
ergies of the states \Ej) and \Ei) relevant to the three 



components as follows; 

V a ; (Ej,Ei) = 

b 



(0,E K + 6E a ), 



rf: (E^Ei) = (^±SE b ,^ T SE b ), 



(18) 
(19) 



v c ; (E„Ei) = (^±SE C ,^ T SE C ). (20) 

The small corrections 5E stem from the small interac- 
tion energy Ed «C Ek- Substituting these eigenenergies 
into Eq. (13), we find that the three components are 
approximated by 
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(21) 
(22) 
(23) 



where we defined u = Ek I '(fcsT). The quantities p a , 
Pb and p c equal to the value of |(-Ej|p x |-E;)| for the case 
of rf , rf and rf , respectively. The explicit form of the 
partition function Z(u) is 



Z{u) = l + 4e-"/ 2 + 4e- 
+ 4e- 4u + 4e- 



Ae- 2u + 8e- 5u/2 
1-'' " +••• , (24) 



which monotonically decreases with rasing u and reaches 
unit for the limit u — > 00. This means that the compo- 
nent rf (u) is a monotonic increase function of u. On the 
other hand, the components rf{u) and n c (u) is convex 
functions giving a maximum at finite u. The conditions 
of u for the maximum of n b and rf are expressed by 



2 



Z'{u) 

Z(u) 

5 Z'(u) 

1 u — u -— 

2 Z(u) 



0, for rf, 
0, for rf. 



(25) 
(26) 



The solutions of the Eqs. (25-26) is estimated as u w 4 
for rf and u w 0.5 for ?7 C . Since the total susceptibility 
X(T) is given by the summation rf +rf '+rj c , it is expected 
that the convex features of rf(u) and rj c (u) cause in the 
bump of the total susceptibility at u w 5 shown in Fig. 
4. 

The argument is clarified by the numerical results 
shown in the inset of Fig. 4, where the u-dcpcndcncc 
of the components for Ed/Ek = 0.1 arc displayed; rf 
(dashed-dotted), rf (dotted), rf (dashed-dotted-dotted), 
together with that of the total susceptibility \ = V a + 
rf + rf (solid) . The component r] b clearly shows a max- 
imum at u w 4, whereas the contribution of rf is negli- 
gible due to the factor e ~ 5u ^ 2 in Eq. (23). As a result, 
the summation rf (u) + rf{u) shows a bump at u = 5.0, 
which is the origin of the anomalous bump of the to- 
tal susceptibility x(T) at the characteristic temperature 
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T* = EkI^b x 0.2. We should note here that, if quan- 
tum rotors are not interacting at all, the component rf 
exactly vanished due to the degeneracy Ej = Ei = Ek / 2 
(See Eq. (19)) and only the component rf is dominant 
for the total susceptibility x(T). This means that the 
total susceptibility is a monotonic function as the same 
as rj a so that the bump does not emerge. The anomalous 
bump of the susceptibility, therefore, manifests the rele- 
vance of the dipolar interaction to the dielectric response 
of quantum rotors. 

V. CONCLUSIONS 

It is important to recall experiments reported in Ref. 
[18], for the dielectric susceptibility of KC1 crystals with 
Li defects. It has been found that the susceptibility 
does not scale linearing with the Li concentration, and 
even becomes smaller with increasing the concentration 
(« 1000 ppm), where the interaction between defects be- 
comes relevant. In addition, a bump of the suscepti- 
bility is observed at about 200 mK for concentrations 
of 200-1000 ppm. These temperature dependences of 
the susceptibility together with the bumps are recovered 
well by our results shown in Fig. 4. Noting that defects 
in both systems move along closed loops and correlated 



each other, it is natural to assume that the similar pic- 
ture holds. For a quantitative discussion, of course, one 
should take into account the effect of potential variation 
hindering the free rotation of Li + , which is caused by the 
Coulomb interaction between a mobile Li + ion and the 
host atoms K + and Cl~ . The problem has been theoret- 
ically investigated in Ref. [19, 20] based on the two-level 
tunneling model. 

In conclusion, we have investigated the quantum dy- 
namics of two dipole rotors coupled via dipolar interac- 
tion. By solving analytically the eigenvalue problem of 
coupled rotors, we have demonstrated the energy spec- 
tra of coupled rotors as a function of dipolar interaction. 
The anomalous temperature dependence of dielectric sus- 
ceptibility is also shown. Our model is so general that 
it should be applicable in a variety of physical context 
relevant to quantum rotors. 
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